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,— ' " The aim of this paper is to illustrate and solve, under a class of generalized Levy models, the pricing problem for a 
^ contingent claim jointly depending on a market asset St and its realized volatility RVt . The contribution of this work 
^ seems to us as being two-fold: 

. . . ■ (%)• It defines a class of asset price models that naturally extends the existing theory of time-changed Levy models; 

(ii). It introduces the class of derivatives paying off on an asset and the volatility in mathematically rigorous terms, 
and provides the necessary machinery for pricing contingent claims of such kind in the presented stochastic 
framework. 

Although these two ideas do represent separate research topics, they will be developed in this article through a 
common discussion. In particular, the elaboration of the first point will be functional to the fulfillment of the second. 



Abstract 

A decoupled time-changed (DTC) Levy process is a generalized time-changed Levy process whose continuous 
and discontinuous parts are allowed to follow separate random time scalings. Disentangling the stochastic time 
change in two components not only provides a powerful unitary framework for the models already present in the 
literature, but in principle opens up new possibilities for the choice of an asset's log-price dynamics. Following 
this new idea, we devise a general martingale structure for the price process, and obtain an inverse-Fourier pricing 
equation for claims paying at maturity a function of the asset value and its realized volatility. Thus with a single 
formula we are able to capture prices for derivatives depending on either an asset or its volatility, but also joint 
payoffs are allowed, like the target volatility option (TVO), which shall be discussed. Numerical computations 
validating our techniques are provided. Notably, the DTC theory allows to incorporate valuation formulae from 
various asset models in a single software implementation. 

Keywords: Derivative pricing; time changes; Levy processes; volatility derivatives; target volatility option; martingale 
representation theory. 
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1.1 Time changes 

The use of Levy models in finance dates back to 1976 when Merton, in [36J, proposed that the log-price dynamics of 
a stock return should follow an exponential Brownian diffusion punctuated by a Poisson arrival process of normally 
distributed jumps. For the first time then the two main shortcomings of the Black-Scholes model, continuity of the 
sample paths and normality of returns, were addressed. Ever since, Levy processes have proved to be a flexible 
and yet mathematically tractable instrument for asset price sampling and modeling. To name just a few, popular 
models are the Variance Gamma (Madan et al., [M]), the CGMY (Carr et al., [7j), and the Normal Inverse Gaussian 
(Barndorff-Nielsen, [3]). 

In finance, the ideas of Levy processes and time changes were linked from the very beginning. Indeed, one of the 
easiest way of producing Levy processes is using the principle of subordination of a Brownian motion W%. If T t is 
an increasing Levy process independent of Wt, then the subordinated process Wr t will still be of Levy type; most 
importantly, its characteristic function is easily recovered from those of T t and X t . Ultimately, this means that asset 
models based on exponentials of subordinated Levy processes can be obtained from the Levy structures of Wt and T t 
with little effort. 

Subordination is the simplest example of a time change. Given a process X t and an increasing process T t , the 
time change of X t according to T t is the process Xt ± ■ Stock price return models depending on time-changed Brownian 
motions have been conjectured for the first time by Clark in PJ. Further theoretical evidence supporting the financial 
use of time-changed models is given by Monroe's celebrated Theorem in |37| . asserting that any semimartingale can 
be viewed as a time change of a Brownian motion. As a consequence, any semimartingale representing the log-price 
process of an asset can be considered as a re-scaled Wiener process. Empirical studies (Ane and Geman, pQ) confirmed 
that normality of returns can be recovered in a new price density based on the quantity and arrival times of orders, 
which justifies the heuristic interpretation of T t as " business time" or " stochastic clock". 

Further advances were later made by Carr and Wu in [8]. They demonstrated that much more general time changes 
are potential candidates asset price modeling. Consequently, they re-interpreted a tremendous amount of models from 
the literature as time changes. Their theory effectively encompasses many of the areas of asset modeling research: 
finite/infinite activity Levy models, stochastic volatility models, and volatility/asset correlations effects (the so-called 
leverage effect). The latter can be captured in the valuation process by introducing a complex measure change inducing 
distributions of a market with no leverage; such a probability measure is known as the leverage-neutral measure. 

However, Carr and Wu's work still leaves out a number of popular possibilities. For example, the widely known 
model by Bates ([2], [3J) and the relatively more sophisticated one by Fang, [T7], are not time-changed Levy as 
understood in [8]. Indeed, in these two models the jump component does not follow the same time scaling as the 
continuous Brownian part. In the Bates model the discontinuities have stationary increments, whereas in the Fang 
model the jump rate is allowed to follow a diffusion of his own. In other words, price models for which the "stochastic 
clock" runs at different paces for the "small" and "big" market movements have already been proposed and tested. 
Statistical analyses (Bates [3J, Fang [TT], Sec. 4) confirm that these models are capable of improved data fitting. 

Motivated by such instances, the natural question arising is then whether it is possible to manufacture consistent 
general time-changed price processes in which the continuous and discontinuous parts of the underlying Levy model 
follow two different stochastic time changes. 

We shall show that the answer is affirmative. Effectively, we extend the idea of [S] of encapsulating different classes 
of models in a single time-changed representation to an even broader range of stochastic processes, and show that these 
are suitable for financial uses. Moreover, we are still able to capture the correlation structures between the "stochastic 
clocks" and the underlying Levy process by a measure change. 

The class of models we propose is obtained by time-modifying the continuous and jump parts of a given Levy 
process X t by two (not necessarily independent) stochastic time scalings T t and Ut satisfying a certain regularity 
condition (Definition 13. 1|) . We call this a decoupled time change. In a formula: 

Xt,u '•= + ^-u t i (l-l) 
where X£ and Xf represent respectively the Brownian and jump components of X t . Since any time-changed Levy 
process Xx t is in particular a decoupled time change, this clearly provides a natural generalization of the theory of [8J. 
How to derive martingale relations for a stock whose log-returns are modeled by Xt,u will be explained in section 2. 
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1.2 Joint asset /volatility derivatives 

Volatility derivatives took off in the late nineties mostly as a market answer to the post-LTCM crisis. Back then hedge 
funds started selling variance swaps at the high market implied volatility levels, and other counterparties were happy 
to buy these products in order to protect their portfolios against possible implied volatility increases. The realized 
volatility of an asset over a period of n days is empirically computed as the standard deviation of the daily log-returns: 



For example, trading a volatility swap maturing in n days necessitates the computation of the fair price for delivery 
of RV n , times a notional amount, in n days. If St is an equity, the short side of an investor's portfolio involving St 
would be hurt by an increase of the terms in the sum (|1.2[) (vega effect), but going long a volatility swap on St may 
offset the corresponding loss. Compared to delta-hedging the position, this strategy has the obvious advantage of 
eliminating the frictions and dangers typically associated with portfolio rebalancing. 

A vast literature on volatility derivatives pricing consequently developed. To name a few, Dupire in |14] was the 
first to note that a contract on the realized variance can be replicated in a model-independent way. Neuberger [38] 
and Demeterfi et al. proposed the replication techniques of a variance swap via the log-contract, an idea which 
was widely expanded in the pricing by robust replication paper [6] by Carr and Lee. For methods relying on specific 
models see for instance Matytsin [33], Elliott et al. [16] . Javaheri [55] et al.. 

A more recent (as of 2008 to our knowledge) market innovation, is that of derivatives and investment strategies 
based on volatility-modified versions of plain vanilla products. Such contracts are able to replicate classic European 
payoffs under a perfect volatility foresight. At the same time, the component of the price due to a vega excess may be 
reduced by using the realized volatility as a normalizing factor. One example of such a product is the target volatility 
option. A target volatility (call) option (TVO) pays at maturity t the value: 



for a strike price K and some other constant a written in the contract. Intuitively, the closer RV t will be to a, the more 
this product will behave like a call option. However, the presence of RVt at the denominator depresses the sensitivity 
of F to a change in volatility. To give just an instance of how this could be beneficial, in times of market hardship a 
TVO allows an investor to take a long position in an option format without having to incur in the high costs caused by 
implied volatility increases. In section 7 we will investigate more the target volatility option and use it as benchmark 
product for numerical validation of our formulae. 

Motivated by the presence of this new class of derivative securities, it is this authors' interest to pursue a compre- 
hensive study, already begun in [42J, of general joint asset/volatility contingent claims. Instruments of this kind can 
be represented, in terms of their underlying stochastic components, in the form F (St, RVt) for a measurable function 
F of two real variables. 

The problem of pricing joint volatility/asset derivatives in general, and the target volatility option in particular, 
has been already addressed by Di Graziano and Torricelli in [T^] for a 0-correlation stochastic volatility model, and by 
the author in [32], for a general stochastic volatility model. We will recover some of the results from the latter in the 
present work. 

1.3 Pricing 

As already stressed, the decoupled time change approach, although having theoretical relevance of its own, has been here 
conceived as pricing-oriented. Traditionally, mathematical finance presents three main research strains for securities 
pricing: numerical simulation, PDE solution approximation methods, and implementation of analytical formulae. 
Regarding the latter, closed versions are rarely available, the Black-Scholes formula being one notable exception. 
However, because of the nature of the distributions for the underlying price densities, analytical representations are 
much more likely to exist for the Fourier transform of the solution (e.g. Heston model, Levy processes). Having 




(1.2) 





a 



(1.3) 
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such formulae at hand, the valuation problem is simply solved by numerically performing the integration inverting the 
transform. 

Fourier methods for contingent claims valuation have been appearing in the literature since the pioneering work 
of Heston, [22 L who was the first to reckon that the inversion of a characteristic solution of a PDE naturally led to 
Option prices. 

Other important Fourier-inversion methods are given in the FFT paper [5] by Carr and Madan, and in Lewis' book 
[32] and successive paper [33]. The advantage of Lewis' approach is that it is valid for all payoffs whose transform is 
known in closed formula, and that it most clearly differentiates between the payoff and model pricing components. 
Furthermore, it provides an easy way for calculating price sensitivities. 

Lewis' original idea of |32| is in fact a no-nonsense approach to claim valuation. Solving a pricing problem is 
equivalent to finding a PDE solution, which is in turn given by the integral of a fundamental solution times the 
terminal condition. In an equivalent dual formulation, the solution can be found by the inversion integral of the 
Fourier transform of the fundamental solution (the " fundamental transform", a characteristic function from the process) 
multiplied by the terminal condition's (payoff) Fourier transform. An analogous argument can be used (see [33 ) when 
the underlying dynamics do not arise from a diffusion problem but are instead directly derived from a Levy model. 

As we will show, it turns out that these ideas have an ad verbatim multivariate extension whenever the payoff and 
the characteristic function from the model include dependence on the quadratic variation of the log-price process. This 
will be the technique we will favour here for contingent claims pricing. 

Obviously, our proofs are in particular valid for plain contingent claims on an asset St] nonetheless, introducing a 
volatility variable in the payoff can be dealt with by a minimal effort. It is also clear that the results we derive equally 
apply for pricing of classic volatility derivatives, like volatility and variance swaps. Indeed, by reverting to the case 
T t = Ut and considering payoff functions in which the first argument does not appear, it is possible to analytically 
price pure volatility derivatives in the standard time-changed Levy processes theory. 

1.4 Theoretical perspective and structure of the paper 

The technical approach to time-changing featured in this paper is rather different from that found in the previous 
research. In our opinion, some of the existing literature does not take into sufficient account the complications that 
may arise in time-changed asset modeling. In a nutshell, these owe to the fact that stochastic time transformations are 
notoriously ill-behaved with respect to the martingale property. Lacking this, it is impossible to even talk about risk- 
neutral valuation. Note that in this regard the fact that X t is a Levy process is irrelevant. For example, a Brownian 
motion re-scaled with the time change given by the first hitting times of the levels t > (an inverse Gaussian process) 
becomes the deterministic process t itself. The key requirement to preserve the martingale property after a time change 
is the uniform integrability of the underlying process X t : a uniformly integrable martingale remains a martingale after a 
time change; this is due to Doob's optional sampling Theorem. When such an assumption is dropped, counterexamples 
are known (see the discussion after Proposition 13. 3p . 

Unfortunately exponential martingales, the fundamental structures of asset modeling, do not normally enjoy the 
uniform integrability property. For instance a geometric Brownian motion, historically the most prominent asset price 
model, is not uniformly integrable on R + . 

Of course, lack of uniform integrability is ultimately not so important in a model where time runs linearly. Indeed, 
for financial purposes we always look at an asset evolution in a finite time period. Still, whenever X t must be subject 
to a time change, we cannot afford to enforce this limitation, otherwise for some finite values of t the set {T t = +00} 
may have positive probability. That is, T t is no longer a valid time change when restricting X t to a bounded time 
interval. 

A shift of paradigm seems advisable. Throughout this paper we hope to provide convincing evidence that (de- 
coupled) time-changed processes are much more conveniently looked at in the bigger picture of the semimartingale 
representation theory. 

The change of approach we propose effectively avoids any loose arguments when time-changing non-uniformly 
integrable martingales. Insisting on an exponential martingale to remain a martingale after a time change is tempting, 
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but simply not possible, because Doob's Theorem may not apply. Instead, we shall identify a one-parameter family of 
exponential local martingales canonically associated with Xx t , whose elements turn out to be true martingales on a 
certain subset 0o of admissible parameters. Determining the shape of Qq is in general a difficult problem. Sufficient 
conditions are in some cases given; examples are the famous Novikov and Kazamaki conditions for exponentials of 
Brownian integrals, which we will recover as particular cases. 

The matter of this article is organized as follows. In section 2 we lay out the assumptions that we shall need; 
in section 3 martingale properties for a decoupled time-changed Levy model will be derived. Section 4 shows the 
fundamental relation between the characteristic function of the joint transition density of (St, (\ogS)t) and the joint 
Laplace transform of (T t ,Ut) as computed in an appropriate measure. Section 5 is dedicated to the derivation of a 
pricing formula for products paying off on St and (log S) t ■ We devote section 6 to explaining some of the models 
present the literature in terms of a decoupled time change, and find the joint characteristic function of (St, (logS')t) 
for each such model. Finally, in section 7 we use our formulae to valuate a TVO under various market conditions and 
asset models, and briefly summarize the work done. The proofs are found in the appendix. 

2 Assumptions and notation 

As customary, our market model is represented by a filtered probability space (ft, T , Ft, P) satisfying the usual con- 
ditions. Throughout the paper we will assume that there exist a money market account process paying a constant 
interest rate r. 

Let St be a non-dividend-paying market asset. St will denote its time-0 discounted value e~ rt St- The total realized 
variance on [0, t] of St is by definition the quadratic variation of the natural logarithm of St, that is : 

TV t := (logS) t = lim V | log S u+1 - log S u \ 2 ■ (2.1) 

The limit runs over the supremum norm of all the possible partitions it of [io , ■ The total realized volatility is the 
square root of the total realized variance. The period realized variance and volatility (or realized variance/volatility 
tout court) are given respectively by RVt = TVt/t and \/RVt (compare with (|1.2|l ). If X t is a semimartingale, by 
taking the limit in (|2.ip it is an easy check that: 

(X) t =X 2 t -2 [ X u -dX u . (2.2) 
Jo 

The scalar product between vectors is throughout indicated by multiplying on the left with the transposed vector 
■ T . If J is an absolutely continuous random variable, we denote by fj(x) its probability density function and by (j>j(z) 
the characteristic function 

<j>j(z) :=E[e- Tj ]. (2.3) 

For a Fourier- integrable function / : C" — > C its Fourier transform will be denoted /. For a complex- valued function 
or a complex plane subset, •* indicates the complex conjugate function or set. 

Unless otherwise stated, when we say that a process is a martingale we mean a martingale with respect to its 
natural filtration. 

The notation for the conditional expectation of a stochastic process X t at time to < t is E to [ • ] . When the 
distribution of a process X t depends on other state variables xt (as in the case of a Markov process) those are 
implicitly understood to be given at time to by xt . If X t is a process admitting conditional laws, the space of the 
integrable functions in the ^-conditional distribution of X t at time to < t is indicated L\ (X t ). Equalities are always 
understood to hold modulo almost-sure equivalence. 

If Xt is an rt-dimensioiial Levy process, the characteristic exponent of X t is the complex-valued function ipx ■ 
C" ->■ C such that: 

E[e l0TXt ] = e*^W (2.4) 
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where ranges in the subset of C" where the left hand side is a finite quantity. 

For a given choice of a truncation function e(x) (that is, a bounded function which is 0(|x|) around 0) the 
characteristic exponent has the unique Levy-Khintchine representation: 

$ x (9) = ifiO - -—- + / (e 10 x -l- i9 T e{x)) V {dx) 1 (2.5) 

where ^ e € MP, S is a non- negative definite n x n matrix with real entries, and v(dx) a Radon measure on R" 
having a density function that is integrable at +oo and 0(|x| 2 ) around 0. We shall normally make the standard 
choice e(x) — a;l| x |<i and drop the dependence of fi on e. (/i, S, v) is then called the characteristic triplet or the Levy 
characteristics of X t . 

A stochastic time change T t is an ^-adapted cadlag process, increasing and almost-surely finite, such that T t is 
an ^-adapted stopping time for each t. A time change of an n-dimensional Levy process X t according to T t is the 
Tr t -adapted process Y t := Xx t - 



3 Martingale relations and asset price dynamics 

The aim of this section is to devise a complex exponential martingale structure naturally associated to a DTC process. 
This construct serves to a twofold purpose. In first place it allows for a DTC-based asset price evolution enjoying 
the martingale property. According to general theory, this in turn enables to postulate the existence of a risk-neutral 
measure that correctly prices the market securities. Secondly, it defines a class of complex measure change martingales 
pivotal to the characteristic function representation of the the next section. 

Let B the vector space of the J-j-supported Brownian motions with drift starting at 0: 

B = {[it + aW t , fi € K, cr > | W t is an ^-adapted Brownian motion and X = 0}, (3.1) 
and J the vector space of the J" t -supported Levy pure jump processes starting at 0: 



J = < cadlag ^-adapted processes X t with stationary indep. increments 



Every Levy process X t starting at can be decomposed as a sum: 



X t = J2(X s -X s -),X = o\. (3.2) 

s<t J 



X t = X c t +X}, (3.3) 

with X£ € B and Xf e J . We shall refer to A t c and Xf respectively as the continuous and discontinuous parts of X t . 

The processes A 4 C and X% will be in the following subject to time changes. Time changes are very general objects, 
so we have to introduce some additional requirements in order for our discussion to proceed towards the construction 
of a martingale suitable for financial modeling. One property we shall assume throughout is continuity with respect to 
the time change. 

Definition 3.1. Let T t be a time change on a filtration Tt- An J^-adpated process X t is said to be T t - continuous 1 if 
it is almost-surely constant on all the sets [T t - , Tt]. 

Obviously, a sufficient condition for T t -continuity is almost-sure continuity of T t . In particular a continuous process 
is T t -continuous if and only if T t is almost-surely continuous. Hence, of particular relevance is the class of the absolutely 
continuous time changes, with respect to which every stochastic process is continuous. Given a pair of instantaneous 
rate of activity processes, that is, two exogenously-given cadlag positive stochastic processes (v t ,u t ), valid time changes 
are given by the pathwise integrals: 

T t = [ v s -ds, (3.4) 
Jo 

U t 



/ u s -ds. (3.5) 
Jo 



6 



The processes vt and ut describe the instantaneous impact of market trading on St at time f, thus providing a sound 
mathematical description of "business activity" over time. 

We are now ready for the main definition. A decoupled time change of a Levy process is the sum of the (ordinary) 
time changes of its continuous and discontinuous parts: 

Definition 3.2. Let X t be a Levy process and T t , Ut two time changes such that T t is continuous and Xf is Ut- 
continuous. Then: 

X T ,u = X c Tt + X* t (3.6) 
is the decoupled time change of X t according to T t and Ut- 

By Jacod |23| . Corollaire 10.12, a first important property of Xt,u is that it is an J-V tA [/ t -semimartingale. To avoid 
degenerate cases, in what follows we always assume T t and Ut to be such that X T and Xfj t are Markov processes 2 . 

We are now ready to define the class of exponential martingales M t (9) canonically associated with Xt,u- The 
following Proposition represents the main theoretical tool of this paper. 

Proposition 3.3. Let X\ be an n- dimensional Brownian Motion with drift and Xf a pure jump Levy process in M. n . 
Let Tj 1 and T 2 be two time changes such that X\ and Xf are respectively and T 2 -continuous. Set X t — X\ + Xf 
and Tt = (T},T?); define Xt. '■= Xi,i + X'L 2 and denote by be the domain of definition of E,\exp(i9 T X?)] . The 
process: 

M t (9,X t ,T t ) = exp {%9 T X Tt - (T^(fi) + T^ x ,{9)) (3.7) 
is a local martingale, and it is a martingale if and only if 9 G Qq, where: 

9 = {9 e 6 such that E[M t (9,X t ,T t )} = 1, Vt > 0}. (3.8) 

For = T t 2 , M t reduces to the process in jHj, Theorem 1, equation (8). It is important to understand that even 
in this simpler case Proposition 13.31 is not a consequence of Doob's optional sampling Theorem (unlike stated in [8], 
Lemma 1) applied to the martingale Z t (9) — exp(i9 T X t — ti/jx(9)), because the latter is not necessarily uniformly 
intcgrable. Indeed, time-transforming a process always preserves the semimartingale property, but the martingale 
property is maintained for uniformly integrable martingales only, as the following important counterexample due to 
C. Sin shows. Choose X t to be a Brownian motion and let dT t = Yf a dt, a > 0, where Y t is the solution of the SDE 
dYt = Y t dX t , Xa = 1. It is proved in [40J that Xx t , i.e. the exponential martingale solving dZ t — Y t a Z t dX t , is a strict 
local martingale 3 . Surprising as it may be, this means that the set Qq may very well trivialize to {0}. That is, some 
choices of time changes may be inherently unsuitable for time-changed asset price modeling. 

Hence, condition (|3.8j) is of critical importance to avoid degeneracies; clearly, Qq can be strictly smaller than the 
whole O, thus care must be exercised when parametrizing objects of the form (|3.7p . price processes in particular (see 
equation (13. 9p further on). Compare also with the interesting examples found in |31j . subsection 3.8. To give a more 
familiar flavour to (|3.8p . observe that in the simple case of Xx t being a one-dimensional Brownian integral, sufficient 
requirements for it to be satisfied are the classic Novikov and Kazamaki conditions ([39], [28j). The set Oq is called 
the natural parameter set. 

Having obtained martingale relations for a stochastic exponential involving Xt,u> the risk- neutral dynamics for a 
DTC Levy-driven asset are given by the obvious choice. Given a one-dimensional Levy process X t and a decoupled 
time change (T t , Ut), we set a price observation at time t > of spot value Sq to be: 

St = S exp(rt + i9 Q X T ,u - T t ^ c x (6 ) - U t ^ d x (9 )) = S e rt M t (9 , X c t + Xf, (U t , T t )) (3.9) 

with 6*o G 6o being such that (|3.9[) is real. By Proposition 13.31 we see that e~ rt St is a martingale, so that St is a 
consistent price process. This is the fundamental asset model we shall use in the rest of the paper. 
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4 Characteristic functions and the leverage-neutral measure 



We shall now study a characteristic function associated with the price process in (|3.9[) . As mentioned in the intro- 
duction, characteristic functions of price processes are an essential mathematical object for contingent claim pricing. 
Compared to the standing literature our framework has two specific aspects. Firstly, the characteristic function we 
shall consider is not that of the (discounted) log-price alone, but it will also incorporate the quadratic variation of the 
log-process. That is, if St is the value of an asset discounted by r, what we are interested in is the joint character- 
istic function of log (St /So) an d TV t . Just as the characteristic function of the log-price allows for the derivation of 
pricing formulae for contingent claims F(St), such a multivariate extension permits valuations of payoffs of the form 
F '(St, TV t ). 

In second place, in our setup we are in presence of a time-changed structure, in both the log-price and (by Lemma 
14.11 further on) its quadratic variation. As pointed out in [5], Theorem 1, there exists a simple expression for the 
characteristic function of a time-changed Levy process Xx t in terms of a Laplace transform of T tl whose argument 
depends on the Levy structure of X t . Unlike the case of simple Levy subordination, such a Laplace transform is not 
calculated in the original distribution of T t , but entails a change of measure which effectively captures the correlation 
between the underlying Levy process and the time rescaling. As it turns out, thanks to Proposition ^. 31 this argument 
makes perfect sense also in the case of a DTC Levy process Xt,u- 

In order to even start the discussion, we must first make sure that the quadratic variation respects the additivity 
and time-changed structure of Xt,u- We have the following "linearity/commutativity property": 

Lemma 4.1. A decoupled time-changed Levy process Xt,u satisfies: 

(X TtU )t = (X c ) Tt + (X d ) Ut = a 2 T t + (X d ) Ut . (4.1) 

That is, the quadratic variation of Xt,u is the the sum of the time changes of the quadratic variations of its continuous 
and discontinuous parts. 

The T t and ^-continuity property of the underlying Levy components plays here an essential role. Without such 
an assumption, this Lemma would simply be false: see Kobayashi, [29J, Example 2.5. 

Since T t and XJ t are of finite variation, the total realized variance of an asset as in (I3.9|) satisfies TVt = —9 2 ) (XT,u)t, 
and as a consequence of Lemma |4. II we have: 

TVt = -9 2 (<7 2 T t + (X d ) Ut ). (4.2) 

Now, assume that St is given by (|3.9p with X t having Levy characteristics (/j,,a 2 ,v). We build a two-parameters 
family of complex-valued change of measure martingales associated with Xt.u hr the following way. Define the 2- 
dimensional processes Ct = (X£,0) and D t — (Xf,i9o(X) d ). Application of Proposition RP1 to Ct + D t guarantees 
that the P-equivalent measure Q(z,u>) << P having Radon-Nikodym derivative: 

z tv 1 

r^-L =Mt((iz8 ,iw6 ),C t +D t ,(Tt,U t )), (4.3) 



is well defined for all z, w <E C such that (iz9o, iw9o) G 8o- By using relation (|4. 2[) and operating the change of measure 
entailed by (|4.3p . the conditional characteristic function of (log(S t /So), TVt) at time t < t is seen to be: 

$t (z, w) =E to [exp(iz log(S t /S ) + iw TV t )\ (4.4) 

=E to [exp(iz(i0 o (^ + X d t ) ~ Tt^xiOo) - U t ^ d x (9 Q )) - iw9 2 (a 2 T t + (X d ) Ut ))} (4.5) 

=E to [eMK^9 ,iw9 ) T (C t + D t ) - Tt(iz4, c x (9 ) + tw9 2 a 2 ) - U t iz^ d x (9 ))] (4.6) 

=EP o [exp(-Tt(0 O M(z - iz) - 2 a 2 (z 2 + iz - 2%w)/2) - U t (iz^ d x (9 ) - ip D (iz9 ,iw9 )))}. (4.7) 

Here E®[ • ]t refers to the conditional expectation taken with respect to the measure Q(z, w). Lastly, we compute 
the characteristic exponents appearing in (|4.7j) . By definition we have: 

V4(#o) - / (e l9oX - 1 - iOoxl\ x \<i)v(dx). (4.8) 
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Regarding ipr>, by the usual properties of the quadratic variation we know that AX t := X t — X t - obeys (AX t ) 2 
A(X) t , and that (X t ) jumps if and only if X t does, whence: 



exp 

izAXf + iw(AX d ) 2 



(4.9) 



ip D (z,w) = logE 
Therefore we calculate further: 

il) D (z, w)= f ( e tzx+lwx2 - 1 - i(zx + W3?)l\*\<x)v{dx). (4.10) 
Jr 

Equivalently, if the Levy density v D (dx x dy) of D t is explicitly known, we have: 

1> D (z, w)= [ (e iz *+ iw y - 1 - i(zx + wy)* KieM < 1 )u D (dx x dy). (4.11) 



The conditional joint characteristic function <$> to {z, w) has been then completely described in terms of the charac- 
teristic triplet of X t = X£ + Xf and of the the joint Q(z, u>)-distribution of T t and Ut- Since we are in a Markovian 
set-up, we can use the more compact notation: 

$ to (z, w) = C® uUt {9 Qi i{z - iz) - 9W(z 2 +iz- 2tw)/2, izi> d x {6 ) - il> D (iz9 , iw9 Q )), (4.12) 

where Cj, Ut indicates the Laplace transform of the conditional joint distribution of T t , Ut at time to with respect to 
Q(z, w), and it is a given function of t — to and initial states T to , Ut ■ As a particular case, if T t and Ut are independent 
$ to factorizes in a continuous and discontinuous part: 

<f> t0 (z,w) = $ c t0 (z,w)<f>i(z,w), (4.13) 

where: 

$ c t0 (z,w) = [exp(-T t (fl /i(« ~ i«) - 2 a 2 {z 2 + iz - 2%w)/2))] (4.14) 

= d^(6 (j,(z - iz) - 9la 2 {z 2 +iz- 2iw)/2), (4.15) 

and 

$l(z,w) =Ef o 2 [exp(-U t {iz^(6 )-il>D(iz9o,iwd Q )))} (4.16) 

= L®l(izi> d x (e )-'tp D (iz9 ,iw9 )). (4.17) 

The measures Qi and Q2 are those induced respectively by the Radon-Nikodym derivatives Mt(iz9o,X£,T t ) and 
M 2 {(iz9 ,iw6 ),D t ,U t ). 

It is of interest emphasizing the meaning of the measure Q(z, w). Consider for a moment the special case of X t being 
independent of T t and Ut- By taking the inner conditional expectation in (|4.4p with respect to <r(T t , Ut) and applying 
the usual properties, one directly obtains ()4.7[) with Q(z, w) = P. Therefore, whenever there is no dependence between 
the time changes and the underlying Levy process, no change of measure is needed in order extract the characteristic 
function of a price process in terms of the rate of activity. Thus, whenever correlation between X t and (T t ,Ut) is 
present, Q(z,w) gives a measurement of how the leverage effect due to this correlation impacts price densities. In 
turn, by the well-known Girsanov's results, this change of measure can be absorbed in the P-dynamics of the asset by 
a suitable alteration of parameters in the distributions of T t and Ut- In accordance with [S], we shall call the measure 
Q(z, w) the leverage-neutral measure. Just as prices in a risky market can be equivalently computed in a risk-neutral 
environment according to a different price distribution, by an appropriate distributional modification valuations in 
presence of leverage effect can be performed in a different economy with no leverage. 

We shall make direct use of the leverage-neutral measure in describing the risk-leverage-neutral joint characteristic 
function $ to (z,w) for some actual models. 
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5 Pricing and price sensitivities 



The characteristic function found in section U serves to the purpose of devising analytical formulae for the valuation of 
European-type derivatives having a sufficiently regular payoff F. In order to do so we will extend the Fourier-inversion 
machinery introduced by Lewis in |33j to our multivariate context. 

Recall that since all the involved processes are Markovian, it makes sense treating $> t(l (z,w) like a Gauss-Green 
integral kernel depending only on some given initial state at to- We have the following 2-dimensional extension of 
Theorem 1 of [33]: 

Theorem 5.1. Let Y t = log St, with St given by 113.9)) . Let F(x,y) € L\ Q (Y tl (Y) t ) for all t > to, be a positive payoff 
function having analytical Fourier transform F{z,w) in a multi-strip Hp = {(z,w) G 0, ct\ < z < a2, Pi < w < 
P2} Q C 2 . Suppose further that <f>t (z,w) is analytical in = {(z,w) G 0, 71 < z < 72, 771 < w < 772} and that 
<&t (z,w) € L 1 (dz x dw). IfSpf] £J, 7^ then for every multi-line: 

Lk lt k 2 = {{x + iki,y + ik 2 ), (x, y) e M 2 } cEf-nSj (5.1) 

we have that the time-to value of the contingent claim F maturing at time t is: 

e -r(t-t ) 

E t0 [e-^-^F(Y t ,{Y) t )} = (5.2) 

iki +00 pik,2-\-oo 

I e-' tw ^ t «S^ z e- r{t - to)lz $ to (-z,-w)F(z,w)dzdw. 

iki—oc Jik 2 —oc 

It is clear that modifying the asset dynamics specifications only acts on $ to , whereas changing the claim to be 
priced influences F only. Also notice that by setting either variable to 0, we are able to extract from (|5.2|) prices for 
both plain vanilla and pure volatility derivatives. In particular, this means that the classic pricing integrals in |32) and 
|33| , are special cases of the equation above whixh can be recovered when F does not depend on the realized volatility 
and when $ to is either obtained from a diffusion Q32J), or a Levy process (|33J). 

In addition, this representation is particularly useful if we are interested in the sensitivities of the claim value with 
respect to the underlying state variables. Let us consider for instance the Delta (sensitivity with respect to the change 
in value of the underlying) and Gamma (sensitivity with respect to the rate of change in value of the underlying) 
of valuations performed through formula (|5.2[) . Call I(r,to, z,w) the integrand on the right hand side of (|5.2[) ; by 
differentiating under integral sign and noting that $ to bears no dependence on St we see that: 

A t := ^E t0 [e-^-^F{Y u (Y) t )} = j^— / / — Z(r, to, z, w)dzdw, (5.3) 

and 



®S ° 47T 2 Jikx-oo Jik 2 ~oo &to 



—r(t—t ) piki+co ^ifc 2 +oo ■ _ 2 



SC2 hi ^ " 1 4_2 II Q2 

uu ^ n J ik\— 00 J ik 2 —oo u t 



Y t := —E to [e- r{t - to) F(Y t ,(Y)t)} = / / —^I( r ,t , z,w)dzdw. (5.4) 



?2 





Mutatis mutandis we can repeat this argument if we want to determine the price sensitivity with respect to the 
quadratic variation (Y) t . Finally, ${ (z, w) could also depend on other state variables (e.g. the instantaneous rate of 
market activity v to ) known at time to- By calling v one such variable we have that: 

V t := ^E t0 [e-^-^F(Y u (Y)t)] = < ~^- (5-5) 
/ / e-' lw ^ t oS- o lz e r{t - to)tz ^-^(-z,-w)F(z,w)dzdw. 

Jiki—oo Jik 2 — ao ®^ 

This is especially well-suited to the case in which <& ta {z, w) is exponentially affine in v, i.e. 

®to(z-, w ) = exp(A(z,w,t - t )v to + B(z,w,t -t )), (5.6) 

when we have: 

-^(-z,-w) = A(-z,-w,t-t )&t (-z,-w). (5.7) 
In the next section we find explicitly $ to for a number of decoupled time-changed models. 
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6 Specific model analysis 



We now describe the underlying DTC Levy dynamics (|3.9[) for some classic asset price processes, and find the corre- 
sponding risk-leverage-neutral characteristic function <I> to (z, u>). In view of the theory developed this far this accom- 
plishes two goals. In first place, it specifies the fundamental piece of information needed to numerically implement the 
results of the previous section. Secondly, it provides convincing evidence that DTCs offer a natural unifying framework 
for a priori different strains of financial asset models (e.g. continuous /jump diffusions, jump diffusions with stochastic 
volatility, Levy processes). 

The models analyzed fall into two main categories: plain Levy processes (Sees. 6.1, 6.3, 6.4) and decoupled time- 
changed finite-activity Levy models 4 (Sees. 6.2, 6.5, 6.6). The first group can be recovered in a straightforward way 
from our framework by setting T t ,Ut — t. The second group retains a non-trivial DTC structure and constitutes an 
ideal testing ground for out theory. 

The approach will be progressing from the more simple model (Black-Scholes), to the more sophisticated one (Fang). 
However, each finite-activity DTC Levy model can be thought of as a particular case of the model in subsection 6.6, 
by considering the appropriate parameters/time change degeneracies. In practice, this means that a single software 
implementation is able to capture every jump diffusion-type model hereby illustrated. Subcases are dealt with by just 
voiding the relevant parameters, as it is clearly seen in the online MATHEMATICA® supplement. 

As one might expect, the value of 9q m Q3.9|) is implicitly given throughout as the natural Oq = —i. The feasibility 
of this choice from a DTC perspective lies in the specific structure of the time changes. Indeed, T t and Ut will always 
be absolutely continuous time changes with rates of activity v t and Ut given by well-behaved diffusions admitting 
unique strong solutions. Essentially, this means that in the cases we cover we always have 0q = in (|3.8[) . 



6.1 Black-Scholes model 

The classic SDE with constant parameters er, r driven by a Brownian motion Wt ■ 

dS t = rStdt + aStdWt (6.1) 

can be trivially recovered from (|3.9p by setting the triplet for the underlying Levy process X t — Xf to be (0, <x, 0) and 
letting T t = t, Ut — 0, so that Xt,u — X t . From ()4.12j) . we have immediately: 

$>t (z,w) = cxp(-(t - t )a 2 (z 2 + iz - 2ot)/2). (6.2) 



6.2 Heston model 

The model proposed by Heston in |22| to explain the leptokurtic feature of price observations is classically represented 
in a risk-neutral measure by the dependent pair of SDEs: 

dS t = rS t dt + ^r t S t dWl ; (6.3) 

dv t = a(0 - v t )dt + ri^dW? . (6.4) 

The parameters a, 9, r\ > need to satisfy the Feller condition of [IS], and the Brownian motions W/, W 2 are correlated 
with (Wf,Wjr) = pt. By a well-known Theorem ([27], Ch. 3, Theorem 4.6), any continuous martingale N t can be 
written as Nt = W/pf t \ for a certain Brownian motion Wt, which implies that the DTC structure of Xx,u m the Heston 
model corresponds to a standard Brownian motion Wt time-changed by: 

T t = v u du. (6.5) 



(i 



Since we are in presence of leverage between X t and T t , to find the characteristic function $ to (z, w) corresponding 
to this model we apply the complex-plane version of Girsanov's Theorem, which gives us the leverage-neutral dynamics 
of vt- This results in changing the dynamics parameters a, 9 of (|6.4p to: 



ipzrj; (6-6) 
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9 Z = a9/a z . (6.7) 

By (|3.9[) this determines <E>j as: 

$ to (z,w) = C\{z 2 /2 + iz/2-iw) =E to L-(- 2 /2+-/2-™)/ t t <^j = £^( z 2 /2 + iz/2 „ iw) (6 8) 

where is given by the dynamics (|6.4p with the parameters (|6.6j) - (|6.7p . The conditional Laplace transform of the 
integrated square root process T t 2 := J Q v*du is well-known 5 analytically (Dufresne |13j). 

The case T t = t reverts back to the Black-Scholes model, when (|6.8p collapses to (|6.2p with cr 2 = i>o- 

6.3 Jump diffusion models 

In their classic works, Merton and Kou ([36J, [30J) proposed to model the log-price dynamics as a finite-activity jump 
diffusion, which in practice has been observed to better fit the smile of option values for short maturity ranges. The 
risk-neutral asset dynamics are given by: 

dS t = rS t - dt + G t S t - dW t + S t - JdN t - n\S t -dt (6.9) 

where Wt is a standard Brownian motion, N t is a Poisson counter of intensity A, and J is the jump size distribution 
J. N t and Wt are assumed to be independent, and the compensator k equals c6j(— i) — 1. In the Merton model 
J is normally distributed J ~ Af(fi,5 2 ), whereas Kou assumed for J an asymmetrically skewed double-exponential 
distribution, that is, the density function fj(x) to be of the form: 

{ e -ax jf x > o 
(6-10) 
Pqe^ x if x < 

for a > 1 , j3 > and p + q = 1 . 

In these models no time change is involved, so Xt^u coincides with the underlying Levy process X t having char- 
acteristic triplet (0, cr 2 , Xfj(x)dx). Thus, to completely characterize Q to (z,w) all we need to do is to determine the 
characteristic exponents and ipn in (|4.8p and (|4.10p . The former is given by the standard theory as: 

V&(0)=A(M0)-1). (6-H) 

In order to compute tfju(z,w) notice that by (|4.9p the process (Xf, (X d ) t ) is just a bivariate compound Poisson 
process of joint density fjji(x,y) and intensity A for each component. Therefore use of (|4.1ip gives: 

ipD(z, w) = A(0j iJ 2 (z, w)-l); (6.12) 

where <f>jj2 (z, w) is the joint characteristic function of J and J 2 . We conclude from (|4.12p that <E> to has the exponential 
Levy structure: 

$ to (z, w) = exp((t - t )(cr 2 (-z 2 /2 - iz/2 + 2iw)/2 + A(0, 7i j2 (z, w) - iz<j>j{-i) + iz - 1)). (6.13) 
Now for the Merton model we have, completing the square: 

<fij tJ 2 (z, w) = / = dx = e 1-2^™ / ^ (6.14) 

J-oo (5V27T J-oo 8\/2-K 

( ip,z-S 2 z 2 /2+in 2 w \ 

exp ^ i-2is 2 w ) 
~ VI - 2iS 2 w ' 
and the integral converges for lm(w) > —l/2a 2 . For the Kou Model we can write: 

<l>j,jz{z>w) = 4>j + ,jI (z,w) + (j>j_ tJ 2_{z,w). (6.15) 
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The characteristic function of the positive part is: 

(j>j +y j2 + (z,w) = e' zx+lwx2 ape- ax dx = ape- L2L ^ 11 I ' e -(v^™^+(«-«)/2v^™) 2 dx ( 6 . 16 ) 

, ,. a / Erfcf-S^V 

=apv 7r e 4i ™ , — 

I 2^vw 

which converges for Im(u>) > 0. The negative part is obtained by a similar calculation. 
6.4 Tempered Levy stable and CGMY 

Another way of obtaining Levy distributions for the asset price is that of directly specifying an infinite activity Levy 
measure v(dx). In this case the asset evolution may be thought of as being solely driven by the infinitely often occurring 
discontinuities of the process, which is therefore going to be of pure jump type plus a deterministic compensator. In 
our language we can represent the underlying stochastics of this type of processes as Xt,u = X t = Xf, with X t being 
a pure jump Levy process of Levy measure v(dx). 

Two instances are the tempered Levy stable process and the CGMY model. Both of these are obtained as an 
exponential smoothing of a stable-type distribution; the second can be viewed as a generalization of the first allowing 
for an asymmetrical skewing between the distribution of positive and negative jumps. The Levy density for a CGMY 
process is: 

Which is well defined for all c+ , c_ , (3 + , 3- > 0, < 2. When a + = a_ one has the tempered stable process. 

For simplicity in what follows we assume a+, a- ^ 0, 1; for such values the involved characteristic functions still exist, 
but lead to particular cases. 

As in subsection 6.3, to fully characterize Q to (z,w) we only need to determine ipx(@) and ipn(z,w). Once this is 
done then (|4.12|) simplifies to yield yet another Levy exponential structure: 

$ to (z,w)=eKp((t-t )($ D {z,w)-izil>%(-i))). (6.18) 

By setting 71 = J'_ 1 xdu(x), the exponent ipx(@) ^ s gi ven ([113; Proposition 4.2) as: 



a + 

p+) 3, j " ' - \y 1 3 j ' 3. 



^{e)= 11+ V{-a + )^c + [[l-^-Y -l + ^+n-a^-cJ I + Pi ' -1-^). (6.19) 



Set 72 = x 2 dv(x); MATHEMATICA® finds the following interesting representation for the positive part ip^ 
of if) D : 

r+°° 2 c . e -P+x 

lp+(z, W ) =izj 1 +iwj 2 + ( e i*x+xwx -\-{i zx + i wx i)) + dx= (6.20) 

Jo x a+ 
22:71 + iw~/ 2 + ic+3 + + [ -w — —2 z — — h iT{-a + ) ) - c+(/3+ - iz) + [ ) (6.21) 



201 2i3+ J ' V w 



i(3+-iz) fl a+ 



w \2 2 



1 — a+ 3 i(/3+ — iz) 



4w 



r(-^],r, 



a+ 1 i(8+ — iz) 



2 ' 2 ' Aw 



Here L is the Euler Gamma and \F\ the confluent hypergeometric function. The strip of convergence of (|6.21j) is the 
set = {(z,w), Im(w) > 0, Im(z) > — 3+}- The determination has a similar expression. 

6.5 Bates model 

The Bates model is a widely known jump diffusion/stochastic volatility model introduced in [2] and [3] to valuate 
foreign exchange options. In the present discussion, it provides us with a first instance of decoupled time change not 
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otherwise obtainable as an ordinary time change. Specifically, the Bates model is a DTC model with a time-changed 
continuous part and a time-homogeneous jump part. 

The dynamics for the asset price and the stochastic volatility are driven by two correlated Brownian motions W} 
and W 2 , and a Poisson process N t of parameter A, combining to the SDEs: 

dS t =rS t -+ y/mSt-dWl + S t - JdN t - n\S t - dt; (6.22) 

dv t = a{9 - v t )dt + ri^tdW?. (6.23) 

The usual choice for the distribution of J is normal 7V(log(l + n) — 5 2 /2, 5 2 ), so that k here represents the average 
percentage jump ratio. The dependence structure and the requirements on the parameters of the driving Brownian and 
Poisson processes are the same as jointly the Heston and Merton models. The Bates model is in effect a combination 
of the two. 

The underlying DTC structure of the Bates model is given by Xx,u = X T + Xf with the characteristic triplet for 
X t being (0, 1, Xfj(x)dx). The time change of the continuous part is, as in subsection 6.2: 

T t = [ v u du. (6.24) 
Jo 

It is straightforward to see that <& ta (z, w) decomposes as: 

$ t0 (z,w) = $ c t0 (z,w)$i(z,w), (6.25) 
where $j o (z,u>) and &f (z, w) are given respectively by (|6.8p and (|6 . 13[) (with a = 0). Thus: 

$ t0 (z,w) =E to L-(^ 2 /2+«/2-i»)/4<d«l ex p(( t _ to )A(^ J j2 ( 0)U ;) -izK- 1) (6.26) 
= C z Tt (z 2 /2 + iz/2 - iw)) exp((i - t a )X(<p JyJ 2 (z, w) - izn - 1). (6.27) 

Expression (|6.27l) is of a completely new kind. Until now we have encountered either exponential Levy mod- 
els, or exponential affine functions arising as solutions of a PDE problem. Here we have a mixture of the two: a 
time-homogeneous jump factor, modeled as a compound Poisson process and a continuous diffusion factor, whose 
characteristic function solves a time-inhomogeneous diffusion problem. Once again we remark the degenerate case 
T t = t, which just yields a Merton jump diffusion with diffusion coefficient -Jvq. 



6.6 Fang Model 

A relatively little known model has been suggested and empirically studied by H. Fang in [TT]. The asset price dynamics 
consist of a jump diffusion with an Heston-type stochastic volatility and a stochastic jump rate, itself modeled as a 
square root process. The risk-neutral dynamics are: 

dS t = rS t - dt + y/viSt- dWl + S t - JdN t - K\ t S t - dt; (6.28) 

dv t = a(6 - v t )dt + rjy/v^dWi; (6.29) 
d\ t = a x (6 x - \ t )dt + vx^dW?. (6.30) 

Both of the diffusion parameter sets must obey Feller's condition. As usual (W t , W 2 ) = pdt. Like in Bates model 
the jumps J are normal with mean log(l + k) — S 2 /2 and variance S 2 /2. Conditional on Xt, N t is a Poisson process of 
jump parameter A t independent of every other process. The Brownian motion Wf is also assumed to be independent 
of all the other random variables 6 . 

This model has a clear DTC Levy structure Xt,u given by T t = v u du, Ut = f X u du, X t having characteristic 
triplet (0, 1, fj(x)dx). In the decomposition (|4.13j) we immediately recognize $£ o as the risk-leverage-neutral charac- 
teristic function of an Heston process of variance Vt, and as that of a compound Poisson process time-changed 
with Ut- The function $t is then given by: 

$ t0 (z,w) = $£ o (z, w)$i (z,w)= C z Tt (z 2 /2 + iz/2 - iw)C Ut {izn - fa^ (z, w) + 1). (6.31) 
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The Laplace transforms of the integrated-square root processes arising from u| and At are known, and the leverage- 
neutral version v* of vt has been given in !6.2l Observe that there is no leverage effect in the jump part because of the 
assumptions on Wf . Finally, notice that the case Ut = t is just the Bates model where the jump activity rate equals 
A . 

7 Numerical testing and final remarks 

7.1 Implementation for the target volatility option 

For validation purposes we numerically implemented the pricing formula (|5.2p on MATHEMATICA® for various 
models, and compared the analytical prices found to a MATLAB® simulation following a naive Euler scheme. 

The target volatility call option is a natural candidate for such testing, being to our knowledge one instance of 
a joint asset/volatility derivatives that has actually been traded. As mentioned in the introduction, a TVO is a 
contingent claim based on underlying asset S paying at the expiration t the amount: 



F{St,RV t ) = 




for a strike price K and a target volatility level a written in the contract. Clearly, the payoff of a TVO will be greater 
than that of the corresponding vanilla call, as long as the volatility realized by the underlying is lower than the target 
volatility. It can be shown (Di Graziano and Torricelli, [12]) that an at-the- money TVO is priced, up to a second 
order error, by the at-the-money Black-Scholes call value of implied volatility a. Thus in some sense a represents the 
volatility prediction an investor attaches to a call option at present time, as opposed to the current market implied 
volatility values. If the latter are too high he may choose to trade in a TVO instead. 
The Fourier transform F(z, w) of a TVO payoff to be used in (|5.2[) is: 

which exists and is analytical in the strip T,p = {(z,w) G C 2 , Im(z) > 1, Im(w) > 0}. 

The model parameters have been taken from the S&P 500 fits of [T7] and are illustrated in table 1. Tables 2 to 
5 summarize the result obtained for four different sets of market and contract conditions. We observe homogeneous 
values for the option prices across the various models for any given data set. Also, an excellent match is exhibited 
with the benchmark Monte Carlo simulation. For other fast and efficient pricing formulae specific to the TVO also 
refer to [T2] , 

Actual computations can be found in the MATHEMATICA® notebook accompanying this paper. A powerful fea- 
ture of the DTC approach can be thereby appreciated. As highlighted in section 6, by voiding the relevant parameters, 
a single coding for equation (|6.3ip is enough to produce prices for several different models. 

7.2 Conclusion 



In this paper we have introduced the concept of decoupled time change, with the intent of accommodating in an 
unitary time change framework models escaping the classic definition of time-changed Levy process. Most notably, 
Bates stochastic volatility /jump diffusion model and Fang's stochastic volatility/stochastic jump rate model are seen 
to have their own, previously unrecognised, time re-scaled Levy structure. 

We obtained martingale relations for stochastic exponentials based on DTC Levy processes, and built up asset 
price dynamics based on such exponential martingales. In order to do so, we had to resort to the semimartingale 
representation of Xt,u as given by its local characteristics. 

Inspired by we proceeded to find a relation between the joint characteristic function of the log-price dynamics 
/quadratic variation and the joint Laplace transform of the time changes. We used this as an integral kernel for pricing 
derivative securities paying off on an asset and/or the volatility it accrues before maturity. In doing so, the idea of 
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leverage-neutral measure and its action on a distribution parameters has been recalled and extended to the present 
set-up. 

Several stochastic models have been analyzed in the DTC framework. In the accounted cases we outlined the 
underlying DTC structure and found the risk-leverage-neutral characteristic function. For numerical comparison and 
validation, we focused on one particular joint asset/volatility derivative, the target volatility option. One important 
contribution from a computational standpoint, is that a single software implementation can output prices for several 
different models. 

8 Appendix: proofs 

The proof of the core Proposition (I3.3[) is based on the theory of martingale representations. In particular, Jacod and 
Shiryaev ([2Ij, Ch. 1, 2) Kuchler and Sorensen ([SI], Ch. 11) and Jacod ([13], Ch. X), provide the necessary tools for 
understanding the key results. As far as this author is aware, this type of approach is relatively unused in the financial 
field, so we briefly illustrate the fundamental concepts. 

Let us introduce the notion of local characteristics of a semimartingale, that play in a general semimartingale 
representation the analogous role of the Levy characteristics. Define the Doleans-Dade exponential of an n-dimensional 
process X t starting at as: 

£{X t ) = e ^-W>/2 J| (1 + AX s ) e - A ^ (8.1) 

s<t 

where X£ denotes the continuous part of X t and the infinite product converges uniformly. This is known to be the 
solution of the SDE dY t = Y t -dX t , Y = 1. 

Next, let e(x) be a truncation function and (at, fit, p(dt x dx)) be a triplet of predictable processes respecting the 
following properties (compare [H], Ch. 2, 2.12-2.14): 

PI. a t has finite variation and values in K"; 

P2. fit is an n x n matrix- valued continuous process of finite variation such that f3 t — f3 s is symmetric positive definite 
almost surely for s < t; 

P3. pidx x dt) is a random measure on R + x R n , almost-surely integrable at +00 and 0(|a;| 2 ) around 0, having the 
property that p({t} x R") < 1; 

P4. Aa t = / e(x)p({t} x dx). 

For 9 G C™, associate with (a t , (3 t , p(dt x dx)) the following complex- valued functional: 

<jr t (0) =i0 T a t -d T p t 0/2 + J (e i0Tx - 1 - i9 T xe(x))p(ds x dx). (8.2) 

This functional is well defined on: 

V = 1 6> G C" such that J J e lf)T x e(x) p(ds x dx) < +00 almost surelyj , (8.3) 

and because of the assumptions above it is also predictable and of finite variation. The local characteristics of a 
semimartingale are then defined as follows: 

Definition 8.1. Let X t be an n-dimensional semimartingale supported on a stochastic basis (f2, J 7 , Ft , P) and let 
T> C C" be like in (|8.3p . Additionally assume that £(^t{&)) 7^ for all 9 G T>. The local characteristics of X t are 
the unique predictable processes (at, fit, p(dx x dt)) satisfying PI, P2, P3, P4 and such that exp(i9 T X t ) /£(^t{9)) is 
a local martingale for all 9 G T>. 

This definition is not the most general, (the condition £(^f t (9)) ^ can be relaxed) but it is the most useful in 
practice. As we will see, the cumulant process ^>t{0) plays the role of a "stochastic compensator" for i9 T X t . The 
measure p is the compensator of the jump measure of X t and it is if and only if X t is continuous. Moreover, X t is 
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a Levy process if and only if its local characteristics are deterministic and of the form (ta, t/3, vdt) for the given Levy 
triplet (a,/3, v), in which case ^t(9) coincides with tipx(9). 

We now explain how the property of continuity with respect to a time change interacts with the local characteristics. 
Indeed, if the condition in Definition 13.11 is met by the involved processes, the local characteristic of a time-changed 
semimartingale are well-behaved, in the sense of the next Theorem. Recall that if B is a Borel space, the time change 
of a random measure p(dt x dx) on the product measure space Q, x B(M. + x R") is the random measure defined by: 

p{dT t x dx)(u,t,B) := p([0,T t (u)] x B)(lo) (8.4) 

for all oj G n, t > and sets B G B(R n ). 

Theorem A. Let X t be a semimartingale having local characteristics (at, fit, p{dx x dt)) and cumulant process ^X (9) 
with domain T>, and let T t be a time change such that X t is T t - continuous. Then the time-changed semimartingale 
Y t = Xr t has local characteristics (q^, f3T t , p(dT t x dx)) and the cumulant process ^>X {®) equals ^^(9), for all 9 G 2?. 

Proof. Kallsen and Shiryaev, Lemma 2.7. □ 
Another important consequence of Definition (I3.1[) is the following: 

Theorem B. Let T t be a time change with respect to some filtration J- t and let X t be a T t -continuous semimartingale. 
For all Ft-predictable integrands H t , we have that Hx t is J- r t -predictable, and: 

H s dX s = ( H T _dX Ts . (8.5) 



This is particularly useful when trying to compute the quadratic variation of a time-changed process. Finally, we 
need a result on linear transformation of a semimartingale and the corresponding change in the local characteristics. 

Theorem C. Let X t be an n-dimensional semimartingale of local characteristics (a t ,fj t ,p(dx x dt)) and let M be a 
linear transformation from R n to R" 1 . Then Y t — MX t is an m-dimensional semimartingale whose local characteristics 
(a Y , f3 Y , p Y (dx x dt)) are given by: 



of = Ma t + / (e(Mx) - Me(x))p(dx x ds) 

JO JS." ,0 g n 

pj = M[3 t M T (8 - 6) 
p Y (B x [0,t]) = p({x e R n \Mx G B} x [0,t]), VB G B(W n ). 

Proof. Eberlein, Papapantoleon and Shiryaev, [TS], Proposition 2.4. □ 
Having at hand these theoretical instruments, we can now proceed to prove Proposition! 



Proof of Proposition \3.!A Let (p, S,0) and (0,0, v) be the Levy triplets of X\ and Xf. Because of the and T t 2 - 
continuity assumption, we can apply Theorem A and we immediately have that the local characterics of X*j X^, 2 are 
respectively (T t V, T/S, 0) and (0, 0, dT?v). 

Now consider the 2n-dimensional martingale Y t — (x^,i,X^, 2 J ■ Clearly, the random measure p Y from the 
characteristics of Y t is supported on the last set of n variables, and its margin on such subspace coincides with dT^v. 
Hence, Theorem C applied to Y t , with M being the n x 2n matrix made by the juxtaposition of two n x n identity 
blocks, implies that the local characteristics of Xt { are T^E, vdT^). 

Now, observe that for every 9 G the exponential £(^ t (9)) is well-defined. Moreover, by [M], Proposition 2.9, 
(i) and |15| . Lemma 2.6, Xr t is quasi-left- continuous 7 , which implies that Vl/^f) has no atoms on the sets {t} x dx. 
Consequently, A^ t (9) = 0; also, since is of finite variation, we have ( x I ,c )f = 0. Hence: 

£(¥*(*))= e*p(¥ t (0)). (8.7) 

In particular this means that £(^t(9)) never vanishes. By using the definition of local characteristics we then see 
that the stochastic exponential M t (9, X t ,T t ) is a local martingale for all 9 G 0. But a positive local martingale 
is a supermartingale, and a supermartingale is a martingale if and only if it has constant expectation. Therefore 
M t (0, X t ,T t ) is a martingale if and only if 9 G ©o- 

□ 
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The proof of Lemma |4. II is a mere application of Theorem B. 



Proof of Lemma \4.1\ Clearly: 



(X T ,u)t = {X c T )t + (X*) t +2(X},X*) t . 



It is clear from the Proposition ^. 3l that Xj, is a continuous process whereas X^j is of pure jump-type. Therefore we 
have (Xj,,Xu)t — 0. Since (X c ) t — (J 2 , the result directly follows from Theorem A, applied to (Xj)t and (X^) t . □ 

We conclude by proving the pricing equation of section 5. The proof by Lewis (|33]. Theorem 3.2, Lemma 3.3 and 
Theorem 3.4) can be easily extended to include dependence of the payoff on (Y) t - 

Proof of Theorem \5.1l By writing the expectation as an inverse-Fourier integral (which can be done by the assumptions 
on F and because <£> to is a characteristic function) and passing the expectation under the integration sign we have: 



E t0 [e- r ^F(Y t ,{Y) t )]=E t0 

p—r{t—to) /"iki +00 pik 2 +oo 



r(t— to) rikx+oa rik^+oo 



4tt 2 



Air 2 



iki — oo J ikv—oo 



ik±—oo J iki — oo 

e -™(Y) t0 5-<» e -'-(t-to)«$ to - w )F(z, w)dzdw. 



(8.9) 
(8.10) 



The only thing to prove is that Fubini's Theorem application is justified. Let N t = M t (8o,X t , (T t , Ut)) be the dis- 
counted, normalized log-price. Define the transition densities p t (x, y) — P(JVt < x, (N t ) < y\ to,Nt , (-^)t )l{x6K,y><A') t0 }j 
and let p t (z,w) be their characteristic functions. For all (z,w) £ ifci,fe 2 we have: 



ik\— oo J iki~oo 



ik\—oo J ik 2 — oo 



p t (—z, —w)F(z, w)dzdw 



«(i0.o S^ z e- r{t - to)lz <S> to (-z, -w) F(z, w)dzdw 



(8.11) 

Pt(—z + iki, —w + ik2)F{z + iki, w + ik2)dzdw. (8.12) 



For x £ R, y > 0, set f{x,y) = e klX k2V F{x,y) and g(x,y) = e klX+k2V p t (x,y). We see that the integrand in 
,12p equals f(z,w)g*(z,w). But now / is L x (dx x dy) because F is Fourier-integrable in (for (z,w) £ 'Ep take 



Re(z) = Re(w) = 0); similarly, g* is L (dz x dw) because of the L 1 assumption on $ to . Therefore, application of 
Parseval's formula (Titchsmarsch, [IT], Theorem 35) yields: 



+oo /•+OQ 



Pt{—z + iki, + ik2)F(z + iki, w + ik^jdzdw = 47r^ 



+oo /•+OQ 



p t (x,y)F(x,y)dxdy 



iTT 2 E t0 [F{N t ,(N) t )} < +oo, 



since F£Ll (N t ,(N) t ). 



(8.13) 
(8.14) 
□ 
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Notes 

1 Jacod [23; uses the terminology Tt- adapted, and Tt- synchronized is sometimes found, but Tt-continuous is also common in the literature, 
and in our view less ambiguous. 

2 In general, time changes of Markov processes are not Markovian. By using Dambis, Dubins and Schwarz's Theorem ([27], Theorem 
4.6) one can build a large class of counterexamples by starting from any continuous martingale which is not a Markov process. 

3 It is well known that in certain circumstances price densities associated through the canonical Feynman-Kac argument with some 
degenerate stochastic volatility models may be strict local martingales. This is a typical phenomenon occurring when the underlying 
diffusion problem does not admit a unique solution, as, for instance, in some GARCH-type models (compare with the counterexample 
provided). In the theory of time changes failure to meet the condition expressed in equation l|3.8|l is precisely the equivalent issue, and the 
case for "volatility explosions" ( |32l Ch. 9) has just recreated in our context. 

Subordinated Brownian models are not covered, since these fail to satisfy the Tt-continuity assumption, which is worth recalling, we 
introduced as a sufficient condition for time-changed asset modeling. In particular, Theorem A does not apply to Brownian subordinatied 
processes, thus in some sense valuation of volatility derivatives under such models represents an ill-posed problem. 

5 In |42| . this author has independently found <E>t in the Heston model by augmenting the SDE system associated to the problem 
l |6.3jl - H6.4|) with the equation dlt = vtdt, and solved the associated Fourier-transformed parabolic equation via the usual Feynman-Kac 
argument. As it must be, the two approaches produce the same answer; therefore 1 15. 21 1 also generalizes equation (3.8) of 1421 . 

6 Although Fang model has this assumption, our general framework notably allows for the presence of a correlation structure between 
Wl, W? and Wf. 

7 See [24], Ch. 1, 2.25, 2.26 
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Table 1: Parameter calibration from the S&P estimations of Fang |17| . Sec. 
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Table 2: 8 to = 100, K = 80, to = 1, * = 2.5, r = 0.06, a = 0.1, 1%, = 0. 





Heston 


Merton 


Bates 


Fang 


Analytical value 


8.5801 


8.5022 


8.5575 


8.8729 


Monte Carlo 


8.5660 


8.5191 


8.5511 


8.8455 



Table 3: S to = 100, K = 120, to = 0.5, t = 4,r = 0.039, a = 0.1, TV to = 0. 





Heston 


Merton 


Bates 


Fang 


Analytical value 


5.0899 


5.7586 


5.1129 


5.0767 


Monte Carlo 


5.2121 


5.8470 


5.0908 


5.0976 



Table 4: S to = 100, K = 100, to = 1.2, t = 1.5, r = 0.072, o = 0.1, TV to = 





Heston 


Merton 


Bates 


Fang 


Analytical value 


1.0155 


1.1038 


1.0258 


1.0532 


Monte Carlo 


1.0119 


1.0962 


1.0259 


1.0560 



Table 5: S to = 100, K = 60, t = 3, t = 6, r = 0.0225, a = 0.1, TV to = 0. 





Heston 


Merton 


Bates 


Fang 


Analytical value 


12.2354 


12.1637 


12.3018 


12.5606 


Monte Carlo 


12.2799 


12.2105 


12.3397 


12.5284 
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